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FOREWORD 


This  renoi't  is  the  second  in  a  series  on  the  mathematical  analysis 
of  transient  thermal  variation  in  clothing.  Thermal  diffusivity  is  the 
physical  parameter  associated  with  transient  heat  conduction  and  is  a 
constant  for  any  homogeneous  material.  However,  a  uniform  during  field 
use  rarely  consists  of  a  single  material,  and  a  method  is  needed  for 
determining  the  heat  conduction  across  a  composite  system  of  non-homo- 
geneous  layers  under  non-equilibrium  conditions. 

This  report  presents  a  mathematical  derivation  of  apparent  thermal 
diffusivity  for  composite  fabric  systems  consisting  of  different  cloth 
layers,  and  demonstrates  a  method  for  experimentally  determining  the 
composite  thermal  diffusivity  based  upon  the  general  equation  of  conduc¬ 
tive  heat  transfer  in  the  unsteady  state. 
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ABSTRACT 


A  mathematical  derivation  of  an  apparent  thermal  diffusivity  for  a 
generalized  composite  fabric  system  consisting  of  different  cloth  layers 
is  presented  in  order  to  provide  a  parameter  for  analysis  of  unidimenslonal 
heat  flow  in  a  clothing  system.  An  experimental  method  is  described  for 
obtaining  the  apparent  thermal  diffusivity  as  a  general  solution  to 
Fourier's  Conduction  Equation,  using  a  boundary  step  change  on  a  semi¬ 
infinite  solid  composed  of  dissimilar  fabric  layers. 


Also  presented  are  several  specific  examples  of  transient  heax.  flow 
in  dual  .layer  systems.  The  apparent  diffusivities  of  each  system  are 
obtained  both  theoretically  and  experimentally,  and  the  results  compared 


graphically. 
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APPARENT  THERMAL  DIFFUSIVITY  OF  COMPOSITE  FABRIC  LAYERS 


1.  Introduction 

In  a  previous  report  (7)  the  concept  of  thermal  diffusivity  was 
applied  to  fabrics  in  order  to  provide  a  parameter  for  the  mathematical 
analysis  of  the  thermal  response  of  cloth  to  non-uniform  boundary  condi- 
a  milcm  during  field  activities  rarely  consist. 

Of  a  single  type  of  material,  and  clothing  for  arctic  use  is  composed 
of  several  dissimilar  clothing  layers.  As  a  result,  the  transient 
response  of  the  clothing  to  thermal  changes  at  the  boundary  is  -  ^not^n 
not  only  of  time  and  temperature  but  the  thermal  properties  of  ^indi¬ 
vidual  layers  as  well.  Therefore,  a  method  is  needed  for  combining 
parameters  of  the  individual  layer  into  similar  parameters  f< or  th» 
composite  layer  in  order  to  mathematically  analyze  transient  heat  trans¬ 
fer  by  conduction  through  the  material, 

Verv  little  can  be  found  in  the  literature  concerning  solutions  to 
transient  heat  transfer  problems  by  conduction  involving  composite  layers. 

In  one  of  his  later  books,  Jakob  (4)  Included  a  small  section  pn+C0®P°®tte 
layers;  however,  he  was  concerned  with  interface  tenperatures_betweenwo 
dissimilar  layers,  and  used  the  methods  of  finite  difference  to  -pprcod- 
mate  his  solution.  Capey  and  McKenzie  (1)  worked  *ith 
in  an  aircraft  wing,  but  this  was  a  finite  composite  slab  which  was  heated 
on  one  side  and  cooled  on  the  other  side  by  convection.  Chen  and 
(2)  reported  on  transient  heat  by  radiation  in  composite  material.  However, 
no  onePhas  reported  a  method  of  combining  parameters  to  measure  unsteady- 
state  heat  conduction  through  a  dissimilar  layer  system. 

The  purpose  of  this  report  is  to  show  that  a  single  parameter  -  the 
thermal  diffusivity  -  is  applicable  to  a  series  of  layers.  The 
tion  of  this  complex  thermal  diffusivity  is  needed  to  permit  a  mathematical 
analysis  of  a  clothing  system  under  transient  thermal  conations. 

2.  Theory 

The  general  equation  which  describes  thermal  conduction  through  a  thick 
slab  of  material  which  is  heated  at  one  boundary  and  cooled  at  the  other  is 

41  (i) 

a*  aj 


where 


T  =  temperature 
t  =  time 

x  =  distance  into  cloth  from  heat  source 


and  is  generally  known  as  the  Fourier  Conduction  Equation  in  One  Dimension 


(3).  The  thermal  diffusivity  c*  is  directly  proportional  to  the  thermal 
conductivity  and  inversely  proportional  to  the  density  and  specific  heat 
of  the  conducting  material.  In  cloth,  <Xis  defined  as  apparent  thermal 
diffusivity  and  was  shown  to  be  a  reproducible  characteristic  of  the  cloth; 
thus  the  apparent  thermal  diffusivity  serves  as  a  parameter  for  adequately 
expressing  the  material’s  transient  heat  conduction  properties. 

In  a  comoosite  layer  of  material,  the  apparent  thermal  diffusivity  of 
the  resultant^will  depend  on  the  thermal  conductivity  k,  and  the  heat  capa- 
XpfoFto.  coabtoK  icata rials.  These  to  turn  must  be  detemtorf  free 
the  individual  conductivities  and  heat  capacities  of  the  various  materVls 
and  then  combined  into  an  *fer  the  composite  layers. 


If  the  composite  cloth  consists  of  n  layers  of  dissimilar  materials, 
the  heat  flow  will  be  assumed  unidirectional,  perpendicular  to  tne  layers, 
as  shown  by  arrows  in  Figure  1.  The  eross-eectional  area  is  t he  same  ,  or 
all  the  materials  and  equal  to  A. 


Assuming  the  following  parameters  for  the  various  layers  are: 


thickness 

fl*  ?2>  ?3#  *  *  *n 

thermal  conductivities 

ki*  k2*  k3*  *  "  hi 

densities 

Pi.  ?2.  63*  ’  * 

specific  heats 

Cj,  c2,  C3.  •  •  cn 

DIRECTION  OF'  HEAT  FLOW 


t 

J 


L 


respectively  for  the  individual  layers 
numbered  from  top  to  bottom. 


a.  Thermal  conductivity  of  the  com* 
poslte  layer 


If  the  temperature  at  the  upper 
edge  of  the  first  material  is  T0,  at 
the  interface  between  the  first  and 
second  material  it  is  Tik  at  second 
interface  T2,  at  third  interface  Tj, 
and  so  on  until  the  lower  boundary 
which  is  Tn,  the  heat  flow  may  be 
considered  the  same  for  each  layer 
and  may  be  represented  by  q.  Then 
for  material  1 


q 


kjA  (Tp  ~  Tj) 
*1 


(2) 


Figure  1.  Cross  sectional  diagram  for  material  2, 

of  ’n«  different  fabric  layers  where  k  A(Tt  -  T->) 

arrows  shew  direction  of  heat  flew,  q  - - - - —  (3) 

T's  are  interface  temperatures  and  *2 

x’s  are  thickness  of  layers. 
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for  material  3» 

kq  A(T2  -  T3) 

4  "  *2 

(4) 

and  so  forth  until  material  n,  where 

kn  A(Tn-l  -  Tn> 

q  = 

(5) 

For  all  the  material 

k  A(T0  -  Tn) 
q  “  xi+x2+x3+* **n 

(6) 

where  k  is  the  thermal  conductivity  of  the  composite  layer, 
equations  2  through  5  for  the  temperature  differences 

q  *i 

To  -  T1  “  k!  A 

Solving 

(7) 

m  m  q  X2 

T1  ‘  T2  a  k^A 

(8) 

s& 

II 

& 

1 

C\J 

(9) 

and  so  forth  until 

m  m  -  q  xn 

An-1  “  An  “  kn  A 

(10) 

Summing  the  n  temperature  differences  inclusive  in  equation  7  through  10 

To  -  Tn  “  A  Gl  +  *1  +  1  *  ‘  *  l)  <U> 

Solving  equation  6  for  the  temperature  difference 

(12) 


(13) 


To 

-Tn 

kA 

(xx 

+  x2 

+  X3  +  • 

*  **n) 

Equating 

equations  11  and  12 

-9- 

kA 

(xi  +  x2  + 

x3 

•  •  < 

-  xn) 

.  a 

A 

\kl 

+  *2 
k2 

+  S+  • 

k3 

.  .3»\ 

kn/ 

Dividing 

both  sides 

by 

a 

A 

8 

*1 

+  2£ 

+  22  +  • 

k 

kl 

k2 

k3 

kn 

where 

=  3 

ci  + 

*2  + 

X3  +  •  • 

’  *n 

(14) 
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In  order  to  simplify  equation  14,  both  sides  can  be  divided  by  E*.  and 
making  the  substitutions: 


*1 

Kl  =  ^5 


r2  m^r 


(15) 


K„  = 


Kn 


where  Ki,  K2  and  so  forth  are  the  proportional  part  conductance  of  the 
individual  layers. 

Substituting  15  in  14 


A  +  A  +  A  +  •  •  •  A 


k  Ki  K2  K3 


(16) 


and 


Simplifying  equation  17 


A  +  A  +  A  +  .  .  •  * 

Ki  K2  K3  Kjj 


n 


T 

E  T 

a*l 


(17) 


(18) 


Equation  18  gives  the  thermal  conductivity  of  the  composite  system  as  a 
summation  of  the  reciprocals  of  the  proportional  thermal  conductances  of 
the  individual  layers. 

b.  The  heat  capacity  of  composite  layers 

In  a  similar  manner  the  heat  capacity  of  the  composite  layer 
system  may  be  obtained  frcs  the  heat  capacities  of  the  individual  layers. 

Multiplying  the  volume  (A  x%)  by  the  density  ,  gives  the  mass  of  the 
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first  layer,  and  multiplying  the  mass  by  the  specific  heat  <*.  gives  the 
heat  stored  in  the  layer  if  its  temperature  is  increased  by  1  degree.  So 
if  the  temperature  is  increased  dT  in  the  time  dt,  the  heat  stored  is 


Ql  -  A  xi  Pi  ci 


dt 


(19) 


For  material  2,  the  heat  stored  when  the  temperature  is  increased  dT  in 
the  time  dt  is 


"-J  =  A  *2  P2  =2  S  (20) 

For  material  3 

Q3  =  A  X3  P3  c3  (2D 

and  so  forth  until  material  n 

Qn  =  A  xnPn  cn  H  (22) 

Then  the  total  amount  of  heat  stored  Q  in  the  composite  layer  is  the  sum 
of  the  heat  stored  in  the  individual  layers.  Then 


dT 


Q  *  V)  °a  =  A(*l  P1  C1  +  *2  P2  c2  +  *3  P3  c3  +’  *  *  *nPn  cn^  it 


(23) 


The  total  heat  stored  Q  in  the  composite  layer  may  also  be  expressed  as  a 
function  of  the  composite  heat  capacity  for  an  increase  of  temperature  dT 
in  the  time  dt  and 


Q  =  A(xi  +  *2  +  x3 


+• 


•X„)P  0  jg 


(24) 


where 


P  =  density  of  the  composite  layer 


and 


c  =  specific  heat  of  the  composite  layer 


Equating  equations  23  and  24 

d.T 

A(xx  +  x2  +.  .  .Xn)pc  Sji  =  A(xiPi  ci  +  X2P2  C2+*  *  ncn  )  Jt 
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reducing 


(Xl  +  x2  +•  •  -Xn)  e°  =  xlPl  C1  +  *2^2  c2*  *  ,xnPn  cn 
and  dividing  through  by  x 


Pc  = 


xlP  1C1  +  *2^2  °2  *  **nPn  °n 


Introducing  a  new  series  of  constants  Sj,  S2,  Sj  .  .  .  Sn 
such  that 


Si  = 

s2  = 

s*  = 


xi Pi  C1 
£* 
X2P2  c2 

E* 

*3?3  c3 

E* 


(26) 


xnP n  cn 


n 


Where  S,  may  be  defined  as  the  proportional  part  of  the  total  heat 
capacity  of  the  composite  system  contained  in  the  first  layer. 

Substituting  26  in  25 

p  c  =  Si  +  S2  +  S3  +•  •  •  Sn 


and 


P°=E  Sa 
a=i 


(27) 


Equation  27  gives  the  heat  capacity  of  the  composite  system  as  a 
function  of  the  proportional  heat  capacities  of  the  individual  layers. 

c.  Thermal  diffualvity  of  the  composite  layer 

The  thermal  parameter  required  as  a  measure  of  transient  temp¬ 
erature  variation  by  conduction  is.<*  as  shown  in  equation  1.  In  cloth, 
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-  *■  - 

SKSMtT-  «- r-«f3  £. 

deriSSfcr1  XoX'uyer.  These  two  parsers  my  now  be  cabined 
into  an  <*  for  a  composite  system  and 

L-r  (28: 

£  ta  £ a 

”,1  SSST  — sr 

specific  heats  of  the  individual  layers. 

In  this  report  only  two  different  cloths  will  be  used  in  each 
experiment  so  that  n  =  2  in  equation  28  and 


cx  = 


K1  K2 


(Kj  +  K2)(S1“+  S2) 


procedure 


The  semi-infinite  solid  method,  as  explained  in  a  previous  report 

w - •~%~o&s^usssi  srsssss* 

ST&'JS:  STS  maintain  it  at  a  constant 

temperature. 

.  ,aT.M  -vlindricsl  mass  of  iron  with  heater  wire  wound  around  the 
A  large  cylindrical  mass  »  The  teraperature ■  of  this 

circumference  was  used  as  a  sic  heat  .  J  T,(Fig.  2)  imbedded 

drum  heater  core  temperature  is^nserf  is 

controlled^th^regulator  D  by  adjusting  the  duty  cycle  of  the  circum¬ 
ferential  heater. 


V  , ,  j; 


Figure  2. 


Wiring  diagram  for  thermal  regulation  of  the  block  and  face 
heaters. 


th.  sisz 

+v«  et*ei^  and  also  minimizes  the  possibility  of  hot  spots  wm.cn 
c^S  J^  if^e  heatS  were  placed*^  the  surface  The  of  t^s 
seccmd  heater,  however,  necessitates  the  inclusion  of  another  heating 
system  with  associated  regulatory  mechanism. 

The  thermistor  T,.  for  the  face  is  set  in  a  hole  in^eco^erplate 

close  to  the  exposed  side.  A  second  thermistor  T2  is +P°S^^tSs  t! 
junction  of  the  steel  block  and  copper  plate.  These  two  the rmi*£ors  Ti 
and  To  are  wired  in  series  and  together  they  maintain  the  face  heater 
labeled  HF  in  the  diagram.  When  the  cold  cloth  is  brought  in  contact 
with  the  copper  plate,  it  causes  the  plate  to  lose  heat,  resulting 
temperature  difference  between  Ti  and  T2.  This  activates  regula  or 
whSh  sends  heat  to  the  heater.  The  influx  of  heat  at  T2  tends  to 
negate  the  transient  effect  at  T«,  thus  reducing  the  possibility  of 
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overshoot  of  the  interface  temperature  due  to  the  finite  time  required 
for  the  diffusion  of  heat  to  the  surface.  Power  is  supplied  to  the  face 
heater  through  variac  V. 


The  heater  and  control  wires  I4  and  I-2  are  brought  up  through  the 
cylinder  and  attached  to  a  terminal  strip  mounted  on  the  tcp  of  the  block. 
The  top  surface  and  the  entire  circumference  of  the  block  are  covered  by 
an  insulating  jacket  of  wood. 


The  thermocouple  used  for  measuring  the  interface  or  boundary  temper¬ 
ature  is  imbedded  into  but  insulated  from  the  exposed  side  of  the  ccpper 
plate.  In  this  manner,  the  thermocouple  makes  thermal  contact  with  the 
copper  as  well  as  the  material  to  be  measured. 


Once  the  step  change  is  realized  at  the  cloth’s  upper  surface,  the 
temperature  distribution  within  the  material  at  some  time  after  t  =  0  is 
dependent  upon  the  thermal  properties,  the  initial  temperature  conditions 
and  the  new  boundary  temperature.  This  relationship  is  expressed  mathe¬ 
matically  by  the  particular  solution  of  equation  1  for  a  semi-infinite 
solid  (3)  which  is 


where 


T  ■  Ts  - 


(30) 


cx-  apparent  thermal  diffusivity  of  composite  layers 
t  =  time 

T0  *  tenperature  of  cloth  at  t  =  o 
Ts  =  temperature  of  the  face 
T  =  temperature  at  distance  x  and  time  t 
£3  =  Gaussian  error  integral  function  of  x  and  t 
x  *  thickness  of  1,2,3  .  •  •  composite  layers 


4.  Re  stilts 

The  first  semi-infinite  solid  used  experimentally  consisted  of  40 
layers  of  cloth:  20  layers  of  nylon  66  alternated  with  20  layers  of  blanket 
cloth.  These  two  materials  were  used  because  all  the  necessary  parameters 
were  available  from  a  previous  report  (7). 

Thermocouples  were  placed  in  the  semi-infinite  solid  down  one  com¬ 
posite  layer  (1  layer  nylon  plus  1  layer  of  blanket  cloth)  and  two 
composite  layers  from  the  top,  giving  x»s  of  ,i?2  and  .344  incbes(tbat 
is,  .437  and  ,874  centimeters  to  keep  the  units  consistent  with  those 
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i„  Table  ft  The  cloth  and  ch^ber  «« 

and  when  U0.4°F.  The  resultant  temperature-tine 

SSJrU’SS-LV-Ul  points  ot  Fijere  3. 


Time  in  minulee 


Fieure  3.  Transient  temperature  by  composite  layers  in  a  semi-infinite 
solid  of  nylon  and  blanket  cloth. 

The  solid  curves  of  Figure  3  were  computed  usingequation  30  and 
T  T  anti  Tts  as  listed  above.  The  value  of  apparent  thermal  diffusiyity 
SiA  was  calculated  from  equation  29  using  the  parameters  of  both  fabrics 
listed  in  Table  1. 

The  computed  parameters  for  the  composite  layer  are.aJfJ 
in  Table  1  to  make  comparison  easy.  The  apparent  thermal  diffus^vity  o 
the  composite  layer  when  converted  to  English  units  is  .0120  in  / 
compared  to  .0225  in2/min  for  blanket  cloth  and  .0030  in  /min  for  nylon 

66. 

Figure  4i»a  plot  of  the  complementary  semi-inf inite  solid.  The 
order  ofthe  fabrics  has  been  reversed  with  a  layer  of 
adjacent  to  the  step  change.  The  purpose  of  this  was  to  find  if  any 
•variations  in  the  transient  beat  flow  occur  because  of  varying  the 
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relative  positions  of  the  cloths.  Agnln  the  circled  dots  are  the  experi¬ 
mental  data  and  the  solid  lines  are  computed  for  an  apparent  thermal 
diffusivity  of  .0120  square  inches  per  minute  and  with  an  initial  ■temp- 
erature  T0  =  45.3 .  a  boundary  temperature  of  Ts  =  102.5*  and  x  aga  n  equax 
to  .172  and  .344  inches. 


Figure  4.  Transient  temperature  by  composite  layers  in  a  semi-infinite 
solid  of  blanket  cloth  and  nylon. 

t 

The  experiment  was  next  tried  with  a  more  practical  relationship: 
two  adjacent  components  of  the  standard  arctic  unifom,  underwear  fabric 
(505t  wool-50*  cotton,  tubular  knit,  natural  finish,  12  oz/yd  )  and 
shirting  (85*  wool-15*  nylon,  0G  108).  ^  properties  of  these  two 
materials  and  their  composite  layer  are  listed  in  Table  II.  The  values 
of  apparent  thermal  diffusivity  in  English  units  are:  shirting  .0123 
in2 /min;  underwear  .0155  in2/ain;  and  the  composite  layer  .0136  in  /min. 


The  temperature-time  curves  for  one  and  two  composite  layers  are 
plotted  in  Figure  5.  The  solid  curve  was  computed  using  equation  30  and 
an  apparent  thermal  diffusivity  of  .0136  in2/mln  with  T0  -  45°F,  Ts  = 
103°F  and  x  =  .145  and  .290  inches.  The  circled  points  are  the  experimen 
t.pi  temperatures  obtained  with  To,  Ts  and  x*s  as  stated  above. 
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TABLE  I 

NECESSARY  CONSTANTS  FOR  COMPUTING  THE  APPARENT  THERMAL 
niFFUSIVITY  OF  A  COMPOSITE  LAYER  OF  BLANKET  CLOTH  AND  NYLON  66 


TABLE  II 

NECESSARY  CONSTANTS  FOR  COMPUTING  THE  APPARENT  THERMAL 
DIFFUSIVITY  OF  A  COMPOSITE  LAYER  OF  UNDERWEAR  FABRIC  AND  NYLON-WOOL  SHIRTING 


CLOTH  .  . 1 

Constants 

Underwear 

Fabric 

x  (cm) 

.204 

.165 

.369 

Proportional  part 

•553 

.447 

1.000 

11  («7S zzps) 

26,58  x  10-5 

33.78  x 

10-5 

14.88  xl0-5 

5  (jcmJJoc  J 

6.14  x  10“2 

4.05  x 

10-2 

10.19  x  10-2 

c<  (cm2/  sec) 

1.33  x  10-3 

1.6?  x 

10-3 

1.46  x  10-3 

of (in2/min) 

.0123 

0155 

.0136 

Figure  5.  Transient  temperature  by  composite  layers  in  a  semi-infinite 
solid  of  underwear  cloth  and  nylon-wool  shirting. 


5.  Discussion 

The  close  agreement  between  the  theoretical  and  experimental  data  of 
Figure  3 ,  4  and  5  shows  that  there  is  an  apparent  thermal  dif fusivity 
associated  with  a  composite  layer  of  dissimilar  materials.  Furthermore, 
this  apparent  thermal  diffusivity  can  be  computed  for  the  composite  layer 
of  cloth  if  the  thermal  conductivities,  thermal  diffusivities  and  thick¬ 
nesses  of  the  individual  layers  are  known. 

The  apparent  thermal  diffusivity  of  the  composite  system  must  be 
known  in  order  to  compute  the  transient  heat  losses  through  a  multi¬ 
layered  uniform,  and  it  can  also  provide  valuable  insight  into  thermal 
problems  in  clothing  design.  As  an  example,  consider  the  two  components 
of  the  arctic  uniform  discussed  above,  the  shirting  and  underwear  fabric. 
The  upper  curve  of  Figure  5  is  for  one  composite  layer  and  shows  the 
thermal  response  with  time  for  an  x  of  .145  inches  (1  layer  shirting  +  1 
layer  underwear)  and  a  composite  apparent  thermal  diffusivity  of  .0135 
in* /min,  These  two  constants  are  combined  to  form  the  upper  limit  to  the 
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integral  of  equation  30.  Inserting  this  upper  limit  into  the  equation 
as  the  exponent  of  e,  it  can  be  seen  that  the  rate  of  transient  response 
it,  determined  by 

-  (27S) 


For  one  composite  layer  of  shirting  and  underwear  fabric, 


2 

(.620) 


An  equivalent  transient  response  offered  by  shirting  alone  would  have  an 
upper  limit  equal  to  -(.620)2  or 


With  an  ©<r equal  to  .0123  in2/min,  the  thickness  of  shirting  required 
to  give  an  equal  transient  response  can  be  obtained  by  solving  for  x,  and 
x  =  .138  inches.  Thus  a  layer  of  shirting  .138  inches  thick  would  have 
an  equivalent  transient  thermal  response  to  a  composite  layer  .145  inches 
thick.  Since  a  comparison  of  insulation  values  under  steady  state  condi¬ 
tions  shows  that  identical  insulation  to  .145  inches  of  composite  layer, 
could  be  equivocated  by  .143  inches  of  shirting  alone,  the  major  portion 
of  the  transient  difference  could  not  result  from  differences  in  "clo 
value".  Although  the  underwear  cloth  may  have  features  which  make  its 
inclusion  in  an  arctic  uniform  quite  desirable  (such  as  the  feel  of  contact 
with  the  skin  and  its  moisture-absorbtive  properties)  a  5>  saving  in  bulk 
and  more  Important,  a  13.7#  savings  in  weight  could  be  realized  by  using 
two  modified  layers  of  shirting  rather  than  a  composite  layer  of  shirting 
and  underwear  cloth. 

The  example  given,  although  not  completely  practical,  does  show  that 
mathematical  analysis  of  transient  thermal  components  could  lead  to  in¬ 
creased  environmental  protection  from  clothing.  By  properly  combining 
materials,  uniforms  can  be  designed  which  have  lighter  weight,  less  bulk 
and  greater  tolerance  time. 
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